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Abstract 

Frumkin-Derjaguin theory of interfacial phase transitions and in particular the concept of 
the disjoining pressure of a planar adsorbed film is reviewed and then discussed in terms 
of statistical mechanical formulations of interfacial phase transitions beyond mean-field. 
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1. Frumkin-Derjaguin Theory of Wetting Phenomena 



The thermodynamics of interfacial phase transitions (IPT) was founded over 70 years ago by 
Frumkin [lj and combined with the concept of the 'disjoining pressure' introduced even earlier 
by Derjaguin in the same context [2] and in the context of adsorption in porous media O H]. 
Frumkin was directly concerned with the stability of a film of liquid adsorbed at a planar 
solid-gas interface, which he treated as the direct generalisation of van der Waals's molecular 
theory of bulk phase transitions to interfacial phase transitions. In particular, he considered 
the thermodynamic stability of a planar fluid film as a function of its thickness by implicitly 
assuming that over-saturated gas can exist in the absence of bulk liquid whenever required to 
enable the system to traverse a van der Waals loop within, for example, an adsorption isotherm. 
The interfacial thermodynamic free-energy is the surface-excess grand potential of the liquid 
film per unit area, "fsG^ historically known as surface tension. When the bulk free energy 
terms are subtracted via appropriate choices of Gibbs dividing surfaces, the conjugate pair of 
thermodynamic variables for isothermal processes (which in bulk are chemical potential p and 
number density p) translate to (field, density) = (p, T), where T denotes the adsorption (excess 
amount of molecules adsorbed per unit area). In the context of wetting films, as considered 
by Frumkin, it is natural to use a suitably defined film thickness £ (rather than T) as the 
effective order parameter, which can be done rigorously by defining the thickness to be the 
adsorption divided by a number density that is fixed for a given isotherm; the obvious choice 
suitable for all but unphysically thick unsaturated films is i = T/Ap where for any particular 
isotherm Ap = p^ — py is the difference between the saturated liquid density and the density 
of coexisting vapour at the chosen temperature T. A transformation from the field p to the 
order parameter £ is then via the usual transformation between the grand ensemble and the 
canonical ensemble: 

rrs G (») = f(i)-»Api, (i) 

where f(£) denotes the interfacial excess Helmholtz free energy per unit area. This statistical- 
thermodynamics immediately yields the second-law of thermodynamics in two equivalent forms: 

djsG(v) = -Ap£dp (2) 
df(£) = pApd£. (3) 

The first of these results ([2]) is just the Gibbs adsorption equation. Frumkin and Derjaguin 
use the mean- field quantity p(£) to define a disjoining pressure 

U(£) = [p - p(£)] Ap , (4) 

where p Q denotes the chemical potential corresponding to an equilibrium film thickness £ Q . 
Note that by construction one can read off the disjoining pressure, for any chosen p Q , directly 
from an experimental adsorption isotherm (i.e. i versus p) and hence we can hereafter write 
p{£) as introduced in Eq. It is also possible to introduce a mean-field potential ?p(£), 
an example of what we now call a Landau potential or in this specific context the interface 
potential, by imagining a constraint sufficient to stabilize planar adsorbed films away from 
equilibrium: 

m^o) = f(£)-p Ap£, (5) 

m - . (a. 

Once a specific p a has been chosen, all values of the film thickness except that for which 
p(£) = p represent the free-energy of metastable or unstable planar films. At p{£) = p Q there 
is always at least one value of £ corresponding to the absolute minimum of tp(£), where the 
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interface potential equals the equilibrium surface tension 7sg(a i o)- 



The Russian School summarise the above statistical-thermodynamics with the following 
transformation of the second law of interfacial thermodynamics: 

y SG = if) + m; dip = -Ud£ ; d 7sG = £dU . (7) 

Here one can note three distinct but equivalent approaches to applying Frumkin-Derjaguin 
theory in practice. Frumkin emphasised plots of the thermodynamic potential versus the 
inverse of the order parameter; i.e. jsg versus l/£ which is the direct equivalent of a p — V 
isotherm in van der Waals' theory of bulk phase transitions. Derjaguin preferred to analyse 
the field as a function of its conjugate order parameter; i.e. fi(£) or equivalently adsorption 
isotherms £(/Jj). In either first-order phase transition between equilibrium films of 

different thickness £ and l\ would appear as a van der Waals loop associated with an equal- 
area Maxwell construction: 

yt-o i / l l \ 

i/ii £ \£o hJ 



to 



d£U(£) = , (9) 



respectively [5]. The final route, based on identifying ip(£), was adopted by physicists as the 
work-horse of density functional theory (DFT) of inhomogeneous fluids [6], although the nature 
and significance of the constraint needed to stabilise a film out of equilibrium has typically 
been treated rather vaguely [7j. In the vicinity of a first-order interfacial phase transition ip(£) 
develops a second minimum directly related to the Maxwell construction for a van der Waals 
loop in U(£): 

i/j(£ ) = Wx) = 7sg(Mo) • (10) 

In the absence of a first-order phase transition in the approach to saturation (// = fJ, S at) every 
chosen value of fi Q < fj, sat corresponds exclusively to an equilibrium system (£ a = £(fx = fJ, ))- 
Then when the above statistical thermodynamics is applied to partial wetting at saturation 
(defining an equilibrium planar film surrounding an adsorbed drop of excess liquid with a 
non-zero value of Young's contact angle) one obtains probably the most well-known expression 
from Frumkin-Derjaguin theory [3] 

d£Il(£) = i/j(£ ) - V(oo) = 7 lv(cos 6 - 1) , (11) 

where the final equality has introduced Young's equation [8] for the contact angle 9. The 
above result is a direct expression of the thermodynamic work needed to quasi-statically thin 
a macroscopically thick film of excess liquid down to its equilibrium thickness. 



Frumkin-Derjaguin theory was used as the thermodynamic basis of the famous analysis of 
Dzyaloshinskii, Lifshitz and Pitaevskii (DLP), who adapted high-temperature quantum field- 
theory to directly calculate the dispersion interaction contribution to the normal component of 
the pressure tensor and hence obtain the excess pressure on the surface due to the presence of 
an adsorbed film. By identifying this quantity with Derjaguin's disjoining pressure DLP were 
able for the first time to show that one could explicitly implement Frumkin-Derjaguin theory 
of interfacial phase transitions within theoretical physics [10] . This work was adopted as 
a foundation stone of modern colloid and interface science, despite the clear warning of DLP 
that their field-theory could only be trusted for macroscopically thick films (films surrounding 
drops with a Young's contact angle of no more than about 0.1 degrees [M])- Much harder 
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to explain is why it was more than 15 years after DLP before the interfacial phase transi- 
tions of Frumkin were 'rediscovered' and named wetting transitions [T2l I13j , Direct use 
of Frumkin-Derjaguin theory constitutes a mean-field theory of interfacial phase transitions 
(MFIPT), because of the constraint required to define a van der Waals loop through values 
of the film thickness where a planar film would be metastable or unstable to long-wavelength 
capillary wave fluctuations (eventually broadening into coexistence between films of different 
thicknesses). Interestingly, the upper critical dimension at and below which ip{£) is renor- 
malised by long-wavelength capillary-wave fluctuations away from planarity is known to be 
d = 3 for short-ranged models [14] and in the presence of dispersion forces naturally occurring 
IPT are believed to be mean-field |15j . 



2. Physical interpretation of the mean-field disjoining pressure 

In his monograph, Frenkel [16] reports that Derjaguin's disjoining pressure was apparently first 
discovered "by D. Talmud and S. Bresler, who noticed in 1931 that two drops of mercury freely 
moving on the bottom of a glass vessel filled with molten paraffin did not coalesce but remained 
separated from each other by a very thin paraffin film", [T7]. Frenkel chooses to review this 
work in the context of the "specific elasticity of thin films". However, the name "disjoining 
pressure" indicates that the Russian School understood that Eq. (j3J) must represent the force 
per unit area exerted by the fluid film on the substrate, minus that which would be present 
for the equilibrium film at /i = fj, Q . In fact, DLP appear to assume that the normal component 
of the pressure tensor is directly proportional to the film chemical potential In which case 
one could envisage direct measurements of the disjoining pressure by placing pressure sensors 
just behind the substrate surface. Mean-field DFT has been used to successfully test this deep 
connection within statistical mechanics [7j. In particular, for a wall-fluid model in which the 
substrate field is represented as a one-body external field v ext {z) the statistical mechanical 
definition of the disjoining pressure must be 

/OQ 
dz[p(z;£)-p(z;Q]v ext (z) , (12) 
-oo 

where p(z;£ ) is the density profile of the adsorbed film at equilibrium (at chemical potential 
p Q ). At saturation (p = fj, sa t) one can replace p(z;£ ) with the wall-liquid density profile 
PSl(z), for the purpose of calculating the disjoining pressure from Eq. ([12]) . DLP express 
Eq. (|12|) in terms of the normal component of a pressure tensor. The full tensor can be defined 
to within an arbitrary path integral (or gauge) from the many-body expression for the total 
momentum J(r,t) [18]. The conservation of momentum (J(r,t)) = leads to a sum rule for 
the equilibrium pressure tensor: 

V^(r) = -p(r)V a v ext (r) . (13) 

In planar symmetry this reduces to a sum rule for the normal component of the pressure tensor, 
as required to justify Eq. (fT2j) : 

PN '(z) = -p(z)v ext \z) . (14) 

The above can be regarded as an expression of Newton's third law; i.e. the force of the fluid on 
the wall is equal and opposite to that of the wall on the fluid. Thus, there are two equivalent 
ways of expressing this force balance, either in terms of the intermolecular forces (known as a 
virial theorem) or, equivalently, in terms of the wall-fluid force as in Eq. (|12p . |19| . 
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Figure 1: Planar mean-field interfacial theory defined from confinement of a liquid film in an 
asymmetric planar slit (see text). 

A potential illustration of the connection between the statistical mechanical disjoining 
pressure and the mean-field film chemical potential is sketched in Fig. 1. Here a liquid film 
is depicted confined between the substrate of interest and a hard-wall boundary. If sum rule 
(|14p is integrated over the hard- wall boundary one finds a conceptually important result: 

Pn(L) = kTp HW (L) , (15) 

where k denotes Boltzmann's constant and the notation highlights the fact that the right-side 
of (fT5"j) only applies at the hard- wall boundary; i.e. p(z = LJ) = phw(L). If even a nm-thin 
film of liquid is adsorbed at a substrate then one must be close to saturation and thus one will 
be able to gently confine the film at a pressure pn(L) close to the vapour pressure. Since the 
second virial coefficient is negative along the liquid-vapour coexistence curve, one must have a 
very low value for the density at the hard-wall (provided of course that the substrate field has 
a range smaller than L). Namely, an incipient drying layer must be present at the hard- wall 
[20] . as sketched in Fig. 1. In this manner the hard- wall is able to suppress capillary- wave 
fluctuations away from planar symmetry, without significantly altering the thermodynamic 
state of the liquid film. In effect, the pore width L is being used to set the thickness t of a 
constrained liquid film. Now consider an infinitesimal change of pore width induced by moving 
one of boundary walls while keeping the other fixed. If the wall that is moved is the hard-wall, 
then standard analysis yields [121 EI] 

3 S\ = -Jd 3 r P (z) dVHW ^- z) 

Pn(L)A , (17) 



8LJ T ^ J ^ dL (16) 



where the integral is over all space (i.e. extends across the hard-wall of area A) and sum rule 
(fT5|) has been inserted to obtain the second line. If the substrate wall was moved instead, then 
a result of identical form to (fT6j) is obtained, except that the hard-wall field vhw is replaced 
by the substrate field (note that the sensible way to carry out this alternative derivation is to 
swap the direction of the z axis so that the substrate is defined to lie at z = L until the final 
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result can be re-expressed back in the geometry of Fig. 1): 




Inserting Eq. (|14p then gives back the previous result (|17p since one need only integrate out 
to the hard- wall or to when the substrate field has fallen to zero (whichever is the least). By 
defining if) = ($7 + p Q V)/A, one can rewrite sum rule (|17p in a form directly analogous to (|6j): 

p N (L)- Po = -^ , (19) 

except that now the disjoining pressure really is defined in terms of the force per unit area 
on the hard-wall; p Q is the value pn{L — > oo). In this interpretation the disjoining pressure 
II(L) is the extra pressure on the film needed to maintain a planar film of width £ ~ L at fixed 
T and fixed p = p Q . Alternatively, the result can be expressed in terms of the conservation 
of momentum normal to the substrate wall; i.e. from (|18p . It is this alternative route that 
is directly expressed in the definition (|12p . To within the validity of replacing L with £, we 
may therefore conclude that the disjoining pressure of Frumkin-Derjaguin theory, defined as a 
chemical potential difference ([6]), is equivalent to the extra pressure on a substrate wall that 
would be induced by an adsorbed planar film of metastable or unstable thickness. This deep and 
confusing correspondence is of great significance theoretically, because expressing the disjoining 
pressure in terms of the conservation of momentum normal to the substrate (|12p allows for 
alternative developments in terms of exact physics l22| 123] . In principal, experiment has 
direct access to the local normal force on the substrate wall and so by expressing inhomogeneous 
fluid phenomena in terms of sum rules involving the exact disjoining pressure, as we shall do 
in the following section, one would be able to directly analyse experimental systems. 



3. Beyond mean-field interfacial phenomena 

The physics discussed so far has either implicitly or explicitly assumed the presence of con- 
straints sufficient to prevent large scale fluctuations away from planar symmetry. Thus, any 
prediction of interfacial phase transitions will be mean-field (MFIPT). To go beyond mean- 
field one must include these fluctuations. There are standard procedures for doing this in 
theoretical physics, which essentially boil down to extending the class of constraints to include 
all physically appropriate classes of symmetry (known as summing over 'all' fluctuations). For 
example, given sufficient numerical resources, one could use DFT to search for solutions of the 
class p(x,z) and p(x,y,z), in addition to p(z). Frumkin-Derjaguin theory would have to be 
similarly extended, since Eq. ([6]) implicitly assumes density profiles of the class p(z) stabilized 
by imposing a non-equilibrium chemical potential p(£ ^ £ Q ). Generalised Frumkin-Derjaguin 
Theory (GFDT) would need to be based on instantaneous disjoining potential profiles II(x) and 
n(x,y), with the sum over fluctuations embodied in formal path-integral definitions defining 
the full statistical mechanical averages [22] 

Wo) = <n> (20) 

£ Q = (!) , (21) 

where the instantaneous film profile is denoted £{x) or £(x, y), depending on the class of broken 
symmetry embodied in the fluctuations. The physical significance of fluctuation physics to 
interfacial phase transitions is heavily dependent on geometry. If the overall geometry is planar 
then three-dimensional systems are not strongly affected by thermal fluctuations, although for 
short-range models (such as Ising models) fluctuations do generate higher-order contributions 
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to interfacial critical phenomena needed to fully explain the observed physics [23]. For quasi 
two-dimensional systems however, even in overall 'planar' geometry, interfacial fluctuations 
are sufficiently dominant to invalidate MFIPT. Let us therefore restrict further discussion to 
the geometry of Fig. 2, while bearing in mind the obvious generalisation to higher-dimensional 
geometry in special cases. Figure 2 depicts a region of the interface at which a molecularly thin 
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Figure 2: A thermal fluctuation of an adsorbed film away from planar symmetry in quasi 
two-dimensional geometry. The instantaneous height of the liquid film along the substrate is 
denoted h(x) = £{x). 

film has fluctuated away from the substrate, known as a 'strong' fluctuation regime. Collective 
fluctuations beyond relevant molecular correlation lengths (typically less than a nm) are much 
slower than molecular time scales and so one expects local mechanical equilibrium. Thus, it 
is no accident that the sketch in Fig. 2 reminds readers of Young's picture of an adsorbed 
drop with a well-defined contact angle arising from the cylindrical (or spherical) cap geometry 
away from the substrate. Laplace's formulation of local mechanical equilibrium demands that 
there be an excess pressure inside the drop required to maintain the non-zero curvature k of 
the film profile. Hence, one can immediately appreciate that a fluctuation of this nature is 
directly linked to an interfacial pressure wave (surface sound) embodied in the x-dependence 
of the instantaneous disjoining pressure TL(x). The lowest level of theoretical physics sufficient 
to describe strong fluctuations requires the replacement of the interface potential ip{£) with an 
interfacial Hamiltonian functional: 

H[£] = 1LV J°° dx | \Jl + {i'f - lj + W[£] . (22) 

The first term in the interfacial Hamiltonian defines the cost of distorting the outer surface of 
the film away from planar symmetry. It might be thought that the use of jlv as the restoring 
force is problematic for nanoscopic fluctuations but one should note that in this description 
the film profile £{x) is already coarse-grained over molecular scales and thus its gradient (the 
dash denotes d/dx) is zero everywhere except where a significant deviation has occurred. The 
final term in (|22p denotes the generalisation of the mean-field interface potential to non-planar 
geometry; i.e. we must use a non-local functional of £(x) instead of simply integrating over a 
mean-field tp(£(x)), as emphasised by Parry and coworkers [23]. W[£] is known as the binding 
potential (functional). 

It is actually possible to use the fundamental statistical mechanics of DFT to regard (|22p 
and hence W[£] as formally exact functionals defined by the unique relationship between the 
field <p = fj, — v ext (x,z) and the instantaneous density profile p(x,z), [6]. The latter can 
be used to self-consistently define the instantaneous coarse-grained profile £{x), as discussed 
below. That is, 4>(x, z) is the formally exact constraint needed to define an arbitrary film 
profile. The final result will then be expressed as a formal sum over all physically realisable 
configurations of the film profile. Note that from this point of view our choice of "fiv as the 
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restoring force for bending the film profile is a formal definition, in the same class as choosing 
a Gibbs dividing surface consistent with macroscopic physics; it effects only the distribution of 
free-energy between interfacial free energy and higher-order line tension contributions (which 
it is not particularly informative to include explicitly in the macroscopic Young-Laplace force- 
balance below). The film profile, for a particular constraint 4>(x,z), is formally given by 
minimization of the functional H[l]; i.e. 8H/81 = 0. The final step is to note that conservation 
of momentum perpendicular to the substrate requires us to identify the instantaneous disjoining 
pressure as the functional derivative 

8£ 

which in any case is immediately recognisable as the required generalisation of Eq. Q in 
Frumkin-Derjaguin theory. Namely, with this identification, setting the functional derivative 
of ([22]) to zero leads directly to the sum rule 



n(aO = -iLV-r 1 = -IlvKx) , (24) 




where k{x) is the coarse-grained interfacial curvature profile. That this realisation of Laplace's 
force balance is a direct consequence of the conservation of momentum normal to the substrate 
can be seen in Fig. 2 if we regard the fluctuation as the formation of a cylindrical-cap section 
of a drop of radius R and in-plane radius r = RsmO, defining Young's contact angle 9. We 
note immediately from (f24"|) that at the top of the drop LT reduces to LT C = jlv/R, the Laplace 
pressure difference across the curved film profile. That is, the pressure under the drop is bigger 
than outside of it by the well-known amount. If sum rule (|24|) is integrated from outside the 
fluctuation to the top of the cylindrical cap then one obtains zero, which is an expression of 
force-balance between the bending of the film profile near the three-phase contact line and the 
'macroscopic' curvature of the cylindrical cap; i.e. 



dx 



ft(x) - ft c = rU c = j LV sin 9 , (25) 



where x c denotes the x-coordinate of the centre of the cylindrical-cap portion. In the macro- 
scopic limit where the three-phase contact line can be treated as the intersection (at angle 9) 
of a planar liquid-vapour interface with a planar substrate then Eq. (|25p reduces to an integral 
over the three-phase contact region alone 

dxfl{x) = 7£V'sin# , (26) 

which is the normal component of Young's force balance through the three-phase contact line 
[231 ESI El ESI EI] ■ Here, the disjoining pressure is balanced by the distortion of the film surface 
as one crosses the three-phase contact line. In colloid and interface science this contribution 
is often equated with the disjoining pressure for planar films (of varying thickness) and is 
separated out from the macroscopic contribution II C ; see e.g. [28]. When the above integral 
is taken over the entire surface, in overall planar geometry (Fig. 2), the implied force balance 
is equivalent to full thermodynamic equilibrium; i.e. (II) = as required for thermodynamic 
equilibrium of a planar film. When the equilibrium geometry is non-planar (e.g. in the pres- 
ence of a macroscopic adsorbed drop) both (|24p and (I26p average to sum rules of identical form 
but now involving (ft) (x) and (k)(x). 
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From the above we can conclude that by construction the interface Hamiltonian theory and 
associated disjoining pressure is fully consistent with the conservation of momentum perpen- 
dicular to the substrate. In field theory, sum rules of this class are known as a Ward identity 
and are equivalent to the condition of translational invariance (of the entire system) in the 
z-direction [22]. It is this physics that enables us to equate the exact disjoining pressure profile 
with the excess force on the substrate surface underneath the film: 

n[^](x) =-J_ oo dz [p(x,z;e)-p(z;Q] , (27) 

where the subscript on the external field emphasises that this is the true substrate field without 
an additional constraint. Note that (j27|) reduces to ()12p when the fluctuations are restricted 
to planar 'breathing modes', in contrast to (|24p . Thus, the general statistical mechanical re- 
sult ([27]) contains more physics than ([22]) . In principle, sensors could be placed just under 
the substrate surface to measure the fluctuating disjoining pressure profile n(x), provided the 
fluctuations were sufficiently slow (to enable a local time average before the fluctuation disap- 
peared). It is thus of some interest to note the intimate connection between a measurement of 
IL(x) and the associated I. First let us define a local contact angle profile (not to be confused 
with Young's contact angle): 

tan0(x) = i'(x) , (28) 
£' 

smO(x) = = . (29) 

'l + {If 



Then the first integral of (|24p allows us to define 9{x) in terms of the disjoining pressure profile: 

9{x) = wT x { — [ X dxfl(x)\ , (30) 

I ILV J-oo ) 



where the local profile is assumed to be planar at the left boundary. Finally, the integral of i 
shows that the local fluctuating profile is fully determined by the disjoining pressure profile 
defining 6{x): 

i(x) = ( dxtanO(x) , (31) 



with the zero of film height defined at the far left boundary. In non-planar geometry (such 
as a system with an equilibrium three-phase contact line) the above results are immediately 
applicable to the final equilibrium profile [23j . For overall planar geometry, the true equilibrium 
film width requires a path integral (|21|) over all fluctuations. In principal either of the profiles 
£{x) or II(x) could be directly measured, but these are not independent measurements because 
they determine each other through the conservation of momentum normal to the substrate. In 
a fully molecular system, this link allows one to uniquely define what is meant by the coarse- 
grained profile £(x), at least in the absence of overhangs [23J. One can go further and define 
correlation function hierarchies from the fluctuating disjoining pressure [22]. It has taken the 
present author his entire working lifetime to grasp the true significance of Derjaguin's disjoining 
pressure to the statistical mechanics of inhomogeneous fluids. 
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